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ABSTRACT 

We report on analysis of 308.3 hrs of high speed photometry targeting the pulsating DA white dwarf 
EC14012-1446. The data were acquired with the Whole Earth Telescope (WET) during the 2008 interna- 
tional observing run XCOV26. The Fourier transform of the light curve contains 19 independent frequencies 
and numerous combination frequencies. The dominant peaks are 1633.907, 1887.404, and 2504.897 fiHz. 
Our analysis of the combination amplitudes reveals that the parent frequencies are consistent with modes of 
spherical degree 1=1. The combination amplitudes also provide m identifications for the largest amplitude 
parent frequencies. Our seismology analysis, which includes 2004-2007 archival data, confirms these identifi- 
cations, provides constraints on additional frequencies, and finds an average period spacing of 41 s. Building 
on this foundation, we present nonlinear fits to high signal-to-noise light curves from the SOAR 4.1m, Mc- 
Donald 2.1m, and KPNO 2m telescopes. The fits indicate a time-averaged convective response timescale of 
t = 99.4 ± 17 s, a temperature exponent N= 85 ±6.2 and an inclination angle of (9,- = 32.9 ±3.2°. We present 
our current empirical map of the convective response timescale across the DA instability strip. 
Subject headings: stars:white dwarfs — stars :oscillations — asteroseismology: general — white dwarfs: 
individual(EC14012-1446) —stars evolution 
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1. INTRODUCTION 

Stellar seismology, also known as asteroseismology , pro- 
vides us with a unique tool for probing the interiors of stars, 
allowing us to study fundamental problems in stellar evolution 
such as energy transport, thermodynamics, and magnetism. 
White dwarfs, the evolutionary endpoint of most stars, are 
particularly important targets for asteroseismology. They are 
structurally simple: an electron degenerate carbon/oxygen 
core surrounded by thin non-degenerate layers of hydrogen 
and helium. DA white dwarfs repr esent ~ 80 percent of all 
white dwarfs (Eisenstein et al. 2006), and they have a nearly 
pure layer of hydrogen on top of a layer of helium. DB white 
dwarfs lack this hydrogen layer, having a layer of nearly pure 
helium overlying a carbon/oxygen core. Lacking substantial 
nuclear reactions, white dwarfs simply cool as they age, pass- 
ing through specific temperature ranges (the DBV and DAV 
instability strips) within which they pulsate. These pulsators 
are otherwise normal objects, so what we learn about their 
structure can be applied to the entire population of stellar rem- 
nants, and further applied to our understanding of their main 
sequence progenitors. 

In this paper, we focus on combining asteroseismology of 
the DAV EC14012-1456 with nonlinear analysis of its light 
curve to provide an empirical description of its convection 
zone. Convection is one of the largest remaining sources of 
theoretical uncertainty in astrophysical modeling. Convection 
is traditionally treate d through mixing length theory (MLT, 
Bohm-Vitense 1958). MLT is a simple, local, tim e indepen- 
dent description first applied to stellar modeling by Biermann 
( 1932). It describes the motions of "average" convective cells 
with a mean size / = aH p , where H p is the local pressure 
scale height and a is an adjustable free parameter. MLT is 
clearly incomplete; we know turbulent flows are complex and 
there is no reason why a should remain constant through- 
out the convection zone of a single star, and certainly not for 
stars of different masses, chemical compositions, or evolu- 
tionary phas es. As an example for wh ite dwarfs, Bergeron 
et al.| ( |1995] > and |Tremblay et al.| ( |20T0) show that model pa- 
rameters such as flux, line profiles, energy distribution, color 
indices, and equivalent widths ar e extremely sens itive to the 
assumed MLT parameterization. Bergeron et al. ( 1995 i find 
systematic uncertainties ranging from 25 percent for effec- 
tive temperatures to 11 percent for mass and radius. The use 
of MLT to treat convective energy transport in white dwarf 
atmospheres thus represents a significant source of physical 
uncertainty in calculations of their atmospheric structure. We 
rely on these models to supply the information about white 
dwarf interiors, masses, and temperatures needed to calibrate 
white dwarf cooling sequences. This in turn produces de tailed 
age estimates for white dwarfs (Ruiz & Bergeron 2001) and 
an estimate of the age of the Galactic disk ( Winget et al .|l987| 
Harris et aL]|2006[ ). An observational test of MLT that leads 
to an improved description of convection is an important goal 
th at will hav e implications beyond the study of white dwarfs. 

Brickhill (1992) was the first to realize that a pulsating 
white dwarf s photospheric flux is modulated relative to the 
flux entering at the bottom of the convection zone by an 
amount that depends on the convection zone's thickness. The 
local convection zone depth is a function of the local effec- 
tive temperature, and this varies during a pulsation cycle. The 
result is a distortion of the observed light curve, with narrow 
peaks and wider valleys (Figure [TJ. Convective light curve 
fitting exploits these nonlinearities to recover the thermal re- 




Figure 1. Portion of a high signal-to-noise SOAR light curve of EC14012- 
1446, showing the narrow peaks and broader valleys indicative of convective 
mixing. The error bars for each point are indicated by the size of the points. 
The dominant pulsation period is RJ 600 s.(l mmi f»l mmag) 

sponse timescale of the convection zo ne. Mathematical de- 
tail s of this technique can be fo und in [Montgomery ( 2005c ) 
and Montg omery et al!] ( |2010b| ). For our purposes, applying 
convective light curve fitting to a target star requires three in- 
gredients: 1) a pulsator with a nonlinear light curve, 2) pre- 
cise knowledge of the star's pulsation frequencies and (/,m) 
values, and 3) high signal-to-noise light curves for use in the 
actual fitting process. 

Asteroseismology provides the tools to identify white dwarf 
pulsation (l,m) values. White dwarfs are g-mode pulsators, 
and each pulsation mode can be described by a spherical har- 
monic of degree I, radial overtone k, and azimuthal number 
m, where m takes integer values between — / and I. Given a 
sufficient sample of excited pulsation modes, we can match 
the observed frequencies with theoretical models. An impor- 
tant diagnostic for g-mode pulsators is the mean period spac- 
ing between modes of the same (l,m) but consecutive radial 
overtone k (e.g., the k and k+ 1 modes). The mean period 
spacing depends mainly on stellar mass. Deviations of indi- 
vidual spacings from this mean value provide information on 
the thickness of the hydrogen and/or helium layers and on the 
chemical profile of the core (Montgomery 2009). We note 
here that the actual value of the radial overtone k can not be 
determined observationally, but must be inferred from theo- 
retical models. 

A second diagnostic is given by the presence of multiplets. 
The multiplet components have the same (k, I) and are further 
described by the azimuthal index m, which takes integer val- 
ues between —I and /. To first order, the frequency difference 
relative to the m = component of the multiplet is given by 
bvkim = -mil (I -Cm), where fl is the rotation period and Q/ 
is a coefficient that depends on the pulsation eigenfunctions 
evaluated in the nonrotating case. In the high-A: asymptotic 
limit for g-modes, Ca ~ !/£(£+ 1), although models predict 
it to vary by w 10 percent over the range of observed periods 
in EC14012. Multiplet structure is a strong indication of a 
mode's / value. We expect a triplet for I = 1, a quintuplet for 
I = 2, and so on. The observed frequency differences (split- 
tings) are a measure of the stellar rotation rate as sampled by a 
given mode. Deviations from equal frequency splitting within 
a single multiplet and changes in splittings from one multiplet 
to the next reveal information about differential rotation and 
magnetic field strength. 

Applying convective light curve fitting to a wide sample 
of pulsating white dwarfs provides an empirical map of how 
the convective response timescale varies as a function of ef- 
fective temperature, and this can be compared with theoreti- 
cal models, both MLT and hydrodynamic. The Whole Earth 
Telescope (WET) is engaged in a long term project to provide 
such a description of convection across the hydrogen atmo- 



Whole Earth Telescope Observations of EC14012-1446 
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Figure 2. The final light curve of EC14012-1446 from XCOV26, representing the combination of all data listed in TablefT] The top five panels cover 400000 s 
(4.63 days) each. The bottom panel displays the entire light curve, with time in days. The dominant pulsation period is ksdOO s. 



sphere DAV (r eff w 11100-12200K) and the helium atmo- 
sphere DBV (7k « 22000-29000K) instability strips. We 
present here our results for the DAV pulsator EC14012-1446. 
Our goals are three-fold: determine accurate frequency and 
(l,m) identifications, obtain several nights of high quality 
light curves, and apply convective light curve fitting to obtain 
EC14012-1446's convective response timescale. In the fol- 
lowing, we present the resulting data set and our analysis of 
the Fourier transforms, discuss the identified frequencies and 
(l,m) values, perform nonlinear light curve fits, and present 
the convective parameters for EC14012-1446. 

2. OBSERVATIONS AND REDUCTIONS 

EC14012-1446 (WD1401-147, B=15.67), is a high ampli- 
tude, m ultiperiodic DAV pulsator dis covered by|Stobie et al.| 
( 1995 1 and observed extensively by Hand ler~et al.| (|2008 ). 
Our XCOV26 observations span March 25 to April 29 2008, 
achieving 80 percent coverage during the central five days 
of the run (Figure [2]). Twenty-seven telescopes distributed 
around the globe participated in XCOV26, fifteen of which 
contributed a total of 71 runs (Table [T) on EC14012-1446. 
The observations were obtained with different CCD photome- 
ters, each with distinct effective bandpasses. We minimize the 
bandpass issues by using CCDs with similar detectors where 
possible and employing a red cut-off filter (BG40 or S8612) to 
normalize wavelength response and reduce extinction effects. 

Standard procedure for a WET run calls for observers to 
transfer raw images and calibration files to WET headquar- 
ters for analysis at the end of eac h night. CCD data reduc- 
tion follows the steps outlined in |Provencal et al.| ( |2009| l. We 
corrected each image for bias and thermal noise, and normal- 
ized by the flat-field. Aperture photometry using the IRAF 
photometry pipeline described by Kanaan et al. (2002| was 
performed on each image, utilizing a range of aperture sizes 



for the target an d selected comparison stars. We used the 
WQED pipeline ( [Thompson & Multeity] 12009} to examine 
each nightly light curve for photometric quality, remove out- 
lying points, divide by a suitable comparison star and cor- 
rect for differential extinction. Our observational technique is 
therefore not sensitive to oscillations longer than a few hours. 
The result is light curves with amplitude variations repre- 
sented as fractional intensity (mmi). The unit is a linear repre- 
sentation of the fractional intensity of modulation (1 mmi «1 
mmag). We present our Fourier transforms (FT) in units of 
modulation amplitude (1 mma = 1 x 10~ 3 ma = 0.1% = 1 ppt). 

The final reduction step combines the individual light 
curves to produce the complete light curve for EC14012- 
1446. We assume EC14012-1446 oscillates about a mean 
light level. This assumption allows us to carefully assess over- 
lapping segments from different telescopes and ident ify and 
correct any vertical offsets. As discussed in detail in|Proven4] 



cal et al. (1 2009 1, we find no significant differences between 
noise le 



the noise levels of amplitude spectra using: 1) the combina- 
tion of all light curves including overlapping segments from 
different telescopes, 2) the combination of light curves where 
we retain only higher signal to noise observations in overlap- 
ping segments and 3) combining all light curves using data 
weighted by telescope aperture. 

The final XCOV26 light curve contains 308.3 h of high 
speed photometry. Our coverage is not complete, and this 
incompleteness produces spectral leakage in the amplitude 
spectrum. To quantify this, we sampled a single sinusoid us- 
ing exactly the same times as the original data. The result- 
ing amplitude spectrum, or "spectral window", is the pattern 
produced in the XCOV26 Fourier Transform by a single fre- 
quency. The Fourier transform (FT) and spectral window of 
the complete light curve are given in Figure [3] 
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* High signal to noise light curves used for light curve fitting. All other observations were 
used to obtain accurate frequency information. 
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Figure 3. Fourier Transform of the 2008 EC14012-1446 photometry observations (note vertical scale in each panel). Solid dots label 14 modes containing 19 
independent frequencies. Arrows indicate a sample of combination frequencies (for example F1+F2). Unlabeled peaks are second order (for example F1+F2) 
and third order (for example F1+F2+F3) combinations. The spectral window is plotted in the last panel. Tables|2]and[3]list exact frequency identifications. 
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Figure 4. FTs of the EC14012-1446 data set subdivided into three chunks of f»185 hr each. The labels in the right of each panel give the time segment covered 
by each chunk. The black points in the bottom panel identify the five largest amplitude modes. Some changes in each FT can be explained as differences in the 
window structure and effective resolution for each chunk. In particular, the final chunk has the least data coverage. We do find evidence of amplitude and/or 
frequency variations during the run. All peaks with frequencies higher than fti 3000 fiHz are combination peaks. 
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FigureS. Spectrogram of the five largest amplitude modes in EC14012- 
1446's FT, arranged in order of decreasing period. The bottom panel gives 
the time-varying spectral window for a sample frequency at 1000 fiHz (the 
pattern given by a single sinusoid sampled at exactly the same times as the 
data). Each panel is normalized to an amplitude of one. We find a slow ampli- 
tude increase in the 1633.907 fiHz (612 s) and 1548.146 /iHz (645 s) peaks. 
In addition, the 2308 /iHz (433 s) peak shows a decrease in amplitude during 
the middle of the run. 



3. FREQUENCY IDENTIFICATION 



3.1. Stability 



Following Proven cal et al.| ( |2009| l, we begin by examin- 
ing the stability of EC14012-1446's pulsation amplitudes and 
frequencies throughout the run. Amplitude and/or frequency 
variations produce artifacts in FTs, and greatly complicate the 
identi fication of intrinsic pulsation frequencies (Kepl er et al.| 
2003 I. We divide the data set into three chunks, each spanning 
«185 h (7.8 d). The FT of each chunk is given in Figure [4| 
Several of the dominant frequencies are consistent to within 
measurement error, but we do find evidence of amplitude vari- 
ations in a number of frequencies. 

We calculated spectrograms for the five largest amplitude 
frequencies to further explore the nature of these amplitude 
variations (Figure BJ. A spectrogram quantifies the behav- 
ior of frequencies and amplitudes as a function of time. Our 
spectrograms are generated by dividing the total light curve 
into multiple 5 day segments, each of which overlap by 4.95 
days. The FT of each segment is a measurement of frequen- 
cies and amplitudes centered on a specific time. Each panel in 
Figure[5]is an amalgam, where each segment FT corresponds 
to a vertical line. The x axis is time, the y axis is frequency, 




4000 
Frequency ((iHz) 



Figure 6. The mean noise as a function of frequency for XCOV26. The data 
set was prewhitened by the 62 largest amplitude frequencies. Our noise is 
frequency dependent, but is near 0.2 mma. This is a conservative estimate, as 
we have probably not prewhitened by all the combination frequencies. 

and amplitude is represented by color and is normalized to an 
amplitude of 1 . The bottom panel in Figure Blgives the pattern 
generated by a single sinusoid sampled with exactly the same 
times as the segment light curves (a time-dependent spectral 
window). The results show a slow 18(7 increase in the ampli- 
tude of the dominant 1633.907 /iHz (612 s) peak, from w 22 
to « 3 1 mma, over the course of the run. We also find a simi- 
lar increase in the 1548.146 /iHz (645 s) peak. The 2308 /iHz 
(433 s) peak shows a 25 percent decrease in amplitude during 
the middle of the run and a large apparent decrease at the end 
of the run. The amplitudes of remaining two frequencies are 
consistent to within 3a, 

3.2. The XCOV26 Fourier Transform 

Armed with amplitude stability information for the five 
largest peaks, we are ready to take a car eful look a t EC14012- 
1446's XCOV26 FT. We use Period04 ( |Lenz & Breger|2005| l 
for Fourier analysis and nonlinear least squares fitting to se- 
lect the statistically significant peaks in the XCOV26 FT. As 
detailed in Provencal et al.| ( [2009) >, we adopt the criterion that 
a peak must have an amplitude at least four times greater than 
the mean noise level in the given frequency range. We de- 
fine "noise" as the frequency-dependent mean amplitude after 
prewhitening by the dominant frequencies. This is a conser- 
vative estimate, as it is impossible to ensure that all of the 
"real" frequencies are removed when calculating the noise 
level. This is certainly true for EC14012-1446, where the 
peaks above w 3000 /iHz are mainly combination frequen- 
cies. Figure [6] displays the mean amplitude, specified as the 
square root or the simple mean power using a boxcar of 100 
/iHz, after prewhitening by 62 frequencies, as a function of 
frequency. Our noise is somewhat frequency dependent, but 
is near 0.2 mma. 

To confirm our uncertainty estimates, we calculated Monte 
Carlo simulations using the routine provided in Period04. 
This routine generates a set of light curves using the original 
times, the fitted frequencies and amplitudes, and added Gaus- 
sian noise. A least squares fit is performed on each simulated 
light curve, with the distribution of fit parameters giving the 
uncertainties. Our Monte Carlo results are consistent with our 
mean amplitude noise estimates of w0.2 mma. 

Our frequency selection procedure involves identifying the 
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Figure 7. Prewhitening of the dominant 1633.907 ^iHz mode in the 2008 FT. 
We begin with the removal of the largest amplitude resolved peak (top panel), 
a careful comparison of the residuals in the next panel with the spectral win- 
dow (last panel), and the subsequent removal of additional frequencies. The 
fourth panel shows the residuals after simultaneously removing 1633.907, 
16233.450, and 1624.017^Hz. The 1633.450 /iHz component (right arrow 
in second panel) is result of amplitude changes. 

largest amplitude resolved peak in the FT, fitting a sinusoid 
with that frequency to the data set, subtracting the fit from 
the light curve, recomputing the FT, examining the residuals, 
and repeating the process until no significant power remains. 
This technique, known as prewhitening, must be employed 
with an abundance of caution, especially since we are aware 
of amplitude and/or frequency modulation in our data set. 
The modulation will create artifacts that masquerade as addi- 
tional frequencies. To illustrate, let us examine the region of 
dominant power at 1633.907 /iHz (Figure [7). Comparison of 
the original FT (top panel) with the spectral window (bottom 
panel) demonstrates that most of the signal is concentrated 
at 1633.907 /iHz. We fit a sinusoid to the data to determine 
frequency, amplitude, and phase, and subtract the result from 
the original light curve. The second panel of Figure [JJ shows 
the prewhitened FT. Careful examination reveals two residual 
peaks (at 1633.450 and 1624.015/iHz identified with arrows) 
that are clearly not components of the spectral window. A 
third peak at 1622 /iHz is part of a window pattern associ- 
ated with 1633.450 /iHz, but the alias pattern is asymmetric. 
Panel 3 shows the results of prewhitening by a simultaneous 
fit of 1633.907 and 1633.450 /iHz. Both frequencies are re- 
moved, and 1622 /iHz is diminished, leaving 1624.015 /iHz. 
We next subtract a simultaneous fit of 1633.907, 1633.450, 
and 1624.015 /iHz, with the results displayed in panel 4 of 
Figure [7] No significant power remains. However, a red 
flag is raised: the separation between 1633.907 and 1633.450 
/tHz is 0.457 /iHz, the inverse of which is the run length. 
The 1633.450 /iHz peak is a manifestation of the amplitude 
changes observed in Figure [5] and is not included in our final 
frequency list. 

Similar analysis of several other frequencies turned out to 
be unexpectedly complex. The power at 1775 /iHz (labeled 
"5" in Figure [3} is unresolved. Figure [8] establishes that 
prewhitening this peak requires 4 closely spaced frequencies, 
with frequency differences very similar to the inverse of the 
run length. This is a clear signature of amplitude and/or fre- 
quency instability. Figure ?? shows the spectrogram of this 
region using the same criteria as the spectrograms in FigureB] 
Unlike the slow amplitude increase observed with the 1633 




1775 1780 1785 1790 1795 

Frequency (u-Hz) 



Figure 8. The 1775 /iHz (563 s) mode in the 2008 FT. This power is un- 
resolved, and requires at least 4 closely space frequencies for prewhitening, 
with frequency differences very similar to the inverse of the run length. This 
is a clear signature of amplitude and/or frequency instability. 

/iHz peak, the 1775 /iHz peak undergoes remarkably sudden 
variations, on timescales of a few days. We also find a a de- 
crease in its frequency of 1.2 /iHz (5 a) over the course of the 
run. A second region of power at 1860 /iHz displays similar 
behavior. Both peaks are labeled in Table [2] 

Our final identifications result from a simultaneous non- 
linear least squares fit of 19 independent frequencies, am- 
plitudes, and phases as well as 68 combination frequencies. 
Combination frequencies are fixed with respect to their par- 
ents but their amplitudes and phases are allowed to vary. Ta- 
ble [2]lists 19 identified independent frequencies, with consid- 
eration to those exhibiting amplitude and/or frequency mod- 
ulation. The list includes a doublet with a splitting of 9.9 
/tHz associated with the 1633.907 /iHz peak (labeled 1 in Fig- 
ure [5J, a triplet at 1887 /iHz with an average splitting of 3.8 
/iHz (labeled 2 in Figure [3]l and a second doublet with a split- 
ting of 3.1 /iHz at 2504 /iHz (labeled 4 in Figure [3). Table|3] 
present s the largest amplitude combination frequencies (see 
Section|4~Tj). 

4. MODE IDENTIFICATION 




Time (d) 

Figure 9. Spectrogram of the 1775 /iHz region in the 2008 FT. This region 
of power exhibits remarkable amplitude changes on a timescale of days. The 
bottom panel gives the time dependent window. The spectrogram amplitudes 
(colors) have been normalized to a value of 1 . 
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Table 2 

Table of XCOV27 Independent Frequencies 



Table 3 

Table of Combination Frequencies 



ID 


Frequency 


Period 


Amplitude S/N Notes 


Frequency 


Period 


Amplitude 


S/N Parents 




(A»Hz) 


(s) 


(mma) 




(MHz) 


(s) 


(mma) 



7 


935.3799 ±0.002 


1069.085 


2.7±0.1 


9 




10 


1021.139±0.002 


979.299 


1.7±0.1 


8 




9 


1104.252 ±0.001 


905.591 


2.2 ±0.2 


10 




1 1 


1155.925 ±0.002 


865.108 


1.9±0.2 


8 




14 


1241.403 ± 0.002 


805.540 


1.2±0.2 


5 




13 


14 18. 369 ±0.002 


705.035 


1.2±0.2 


5 




3a 


1521.575 ±0.002 


657.214 


2.2 ±0.2 


9 




3 


1548. 146 ±0.001 


645.933 


7.9±0.1 


32 




la 


1624.015 ±0.003 


615.758 


3.1±0.2 


13 




1 


1633.907 ±0.001 


612.030 


25.7 ±0.1 


104 




5 


1774.989 ±0.100 


563.384 


7.2±0.3 


9 




8 


1860.248 ±0.400 


537.563 


6.4±0.4 


11 




2a 


1883.555 ±0.003 


530.911 


1.5 ±0.2 


5 




2 


1887.404 ±0.001 


529.828 


20.7 ±0.1 


79 




2b 


1891.142±0.002 


528.781 


3.8±0.2 


14 




6 


2304.745 ±0.001 


433.887 


4.7 ±0.1 


23 




4 


2504.897 ±0.001 


399.218 


12.7±0.1 


71 




4a 


2508.060 ±0.002 


398.715 


2.1 ±0.2 


12 




12 


2856.155 ±0.002 


350.121 


2.0 ±0.2 


9 





Frequencies displaying long timescale amplitude modulation 
** Frequency IDs with small letters indicate components of the 
same multiplet. 

Our ultimate goal is to use EC14012-1446's nonlinear light 
curve to extract the star's convective parameters. The require- 
ments for convective light curve fitting include precise knowl- 
edge of the target star's frequencies and amplitudes, as well as 
a good idea of the (/,m) values for these frequencies. The pa- 
rameter space to be explored to determine the best convective 
light curve fit is greatly constrained by knowledge of (l,m) 
indices for the excited pulsations. Montgomery's first appli- 
cation of thi s technique fo cused on the DA pulsators G29-38 
and GD154 (Montgomery 2005a). In the data sets used, these 
objects were dominated by a single pulsation mode, so the 
number of potential (l,m) identifications was small e nough 
that all possibilities could be explored. Montgomery et aL] 



(2010b) recently expanded this work to the multiperiodic pul- 
sator GD358, a well-studied object with detailed (l,m) identi- 



fications of i ts puls a tion modes (Provenc al et al.| ( [2009[ >; [Met^] 
calfe et al.| p000l >; |Winget et al | ( |l9^4| >). WMfeXCOV26 
provided precise frequencies and amplitudes for the modes 
in EC14012-1446 (see Table [2}, we lack prior (l,m) identi- 
fications. Even limiting ourselves to /=1 and 2, as these are 
the spherical degrees most often observed in pulsating white 
dwarfs, this yields a total number of possible combinations of 
order (2/ + 1) 19 , a very large number indeed. In this section, 
we will rely on analysis of the combination frequencies and 
the support of asteroseismology to constrain (l,m) identifica- 
tions in EC 140 12- 1446. 

4.1. Combination Frequencies 

Combination frequencies are typically obs erved in the FTs 
of mo derate to large amplitude pulsators (e.g., Provencal et al. 
[20091 IDolez et alj2006[ |Handler et al.|2002[ >; they are iden- 
tified by their relationships, which must be exact within mea- 
surement error. Combination frequencies can be integer mul- 
tiples of a single parent (harmonics), or sums (or differences) 
of any two modes. These frequencies are not independent, but 
result from nonlinear effects most likely associa ted with the 
surface convection zo ne (Brickhill 1992, Brassard et al.|1995j 
|Wu|2001||lsing & Koester|2001| ). |Wu||2001| i shows that the 



253.466 ± 
2569.286 ± 
2738. 153± 
3096.255 ± 
3182.045 ± 
3267.794 ± 
' 3408.939 ± 
: 3409.912± 
3435.538 ± 
: 3494.064 ± 
3521. 306 ± 
3525.050 ± 
3774.81()± 
3938.646 ± 
4026.445 ± 
4053.043 ± 
4138. 800± 
4192.150± 
: 4279.907 ± 
' 4365.040 ± 
: 4365.608 ± 
4392.348 ± 
4396.055 ± 
4743.561 ± 
4809.633 ± 
5009.696 ± 
5155. 187± 
5361.052± 
5408.698 ± 
5913. 800± 



0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.001 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 



3945.302 
350.121 
365.210 
322.971 
314.263 
306.015 
293.346 
293.263 
291.075 
286.152 
283.986 
283.684 
264.914 
253.894 
248.358 
246.728 
241.616 
238.541 
233.650 
234.464 
229.063 
227.669 
227.477 
210.812 
207.916 
199.613 
193.979 
186.531 
184.887 
169.096 



2.9 ±0.2 
0.6±0.2 
0.9±0.2 
0.6±0.2 
2.7 ±0.2 

4.3 ±0.2 
1.6±0.2 
1.7±0.2 

2.0 ±0.2 
0.8±0.2 
7.0±0.2 
0.9 ±0.2 

2.1 ±0.2 
0.9±0.2 
1.4±0.2 
0.7 ±0.2 
2.0±0.2 
0.7 ±0.2 
2.6 ±0.2 
0.9 ±0.2 
1.3±0.2 
1.5±0.2 
1.8±0.2 
0.9±0.2 

2.2 ±0.2 
0.7 ±0.2 
2.2 ±0.2 
1.2±0.2 
1.3±0.2 

1.4 ±0.2 



9 
4 
4 
3 

12 

13 
5 
4 
6 
4 

33 
5 

12 
5 
7 
4 
8 
4 

10 
4 
5 
7 
8 
5 

12 
4 
8 
6 
6 
7 



2-1 
1+7 
1+9 
2f3 
1+3 
2fl 
6+5 
1+5 
2+3 
1+8 
1+2 
l+2b 
2f2 
1+6 
4+3a 
3+4 
1+4 
2+6 
4+5 
4+8 
4+8 
2+4 
2b+4 
2+12 
4+6 
2f4 
1+1+2 
3+10 
3+3+1 
1+4+5 



* Combination including parent with amplitude modulation 

observed amplitudes of combination frequencies depend on 
geometric factors such as the (l,m) indices of the parent(s) 
and the inclination of the pulsation axis to the line of sight. 
EC14012-1446's FT contains a rich distribution of combina- 
tions that involve 13 of the 19 independent frequencies listed 
in Table [2] Our goal is to exploit the geometric sensitivities to 
provide (l,m) constraint s for these frequ encies. 

|Wu| pOUT) and |Yeates et al] ((2005 ) lay the foundations 
for our analysis by describing analytical expressions for the 
predicted amplitudes and phases of combination frequencies. 
These quantities depend on the inclination angle 8 of the pul- 
sation axis to the line of sight, the (/,m) indices of the parent 
mode(s), the amplitudes of the parent mode(s), and param- 
eters describing the convection zone. To minimize the de- 
pendence on the convective parameters and focu s on the geo- 
metric factors (l,m) and 6, our analysis follows Yeates et al. 
(2005) and considers only combination frequencies that are 
the sum of two parent freque ncies . We combine a genetic al- 
gorithm ( Charbonneau| 1 995| l with |Wu ( 2001 l's formulae, and 
incorporate an improved treatment of limb darkening ta ken di - 
rectly from the models of Koester (Montgomery et al. 2010b). 
For a single run of the code, the best simultaneous fit to the ob- 
served amplitudes of the parent and combination frequencies 
utilizes multiple generations and minimizes the root-mean- 
squared residuals, Res rms , between the predicted and observed 
amplitudes. In practice, we run the code 1000 times, and se- 
lect the solutions having Res lms below a limiting value. This 
process produces a sample of best-fit solutions whose distri- 
bution provides information on the range of values allowed 
for these parameters. 

To test that this approach recovers know n input, we used 
the nonlinear light curve fitting code of Montgomery et al. 
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Figure 10. Probability distribution (from 0-1) of /, in, 8, and A for the 
1633.907 fiHz (612 s) and the 2504.897 /iHz (399 s) variations. The solu- 
tions are from individual runs with Res lms < 0.38 mmi. The top panel gives 
the probability distributions of 8 and A„ found for all frequencies, with pre- 
ferred values of 8 ~ 32° and A ~ 20, respectively. The middle and lower 
panels give the / and m identification distributions for the 1633.907 /tHz (612 
s) and 2504.897 /iHz (399 s) mode. The amplitudes of EC14012-1446's ob- 
served combination frequencies argue that 1633.907 /tHz is 2=1, m=0, while 
2504.897 /tHz is Z = l,m=l. 



(2010b ) to generate a synthetic light curve based on EC 140 12- 
446's parent frequencies. The nonlinear light curve fitting 
code is discussed in more detail in Section[5] For our purposes 
in this test, the importance of using synthetic light curves gen- 
erated by this code rather than a simple simulation employing 
multiple sine functions is that the synthetic light curve will in- 
clude combination frequencies produced via nonlinear effects 
due to convection (Montgomery et al. 2010b). Our purpose is 
to recover the input frequency identifications using these com- 
bination frequencies. We assigned reasonable (l,m) values to 
the input frequencies, chose values for 8, the time-averaged 
convective response time To, and A , where A describes the 
response of the stellar material to the pul sations and includes 
a bolometric correction factor (Wu 2001 ), and added noise. 
Using the method outlined in the previous paragraph, we suc- 
cessfully recovered (/,m) for all large amplitude input parent 
frequencies, while experiencing some disagreement with low 
amplitude parents. This is because low amplitude modes have 
even lower amplitude combination frequencies, which are dif- 
ficult to detect and so are not as numerous. Since our method 
treats both the high amplitude parent modes and the low am- 
plitude combination frequencies equally in the fits, this has 
the effect of de-emphasizing the importance of low amplitude 
combination frequencies. 
Figure 10 shows the resulting probability distribution of I, 



m, 9, and A Q values for fits having Res lms < 0.38 mmi, for 
EC14012-1446's 1633.907 /iHz (612 s) and 2504.897 /iHz 
(399 s) variations. For 1633.907 /iHz (612 s), analysis of 
combination amplitudes strongly argues that the variations are 
best represented as spherical degree 1=1 and azimuthal index 
m=0, with an inclination angle of 30 ± 10°, and a value of 
A = 20 ± 2. We find nearly identical results for 1 887.404 /iHz 
(529 s) and 1548.146 /iHz(646 s). For 2504.897 /xHz(399 s), 
we find the same distributions for 8 and A a , but in this case 
(/,m)=(l,l) is strongly preferred. As we found previously, 
the statistical significance of (l,m) identifications determined 
by combination amplitudes is amplitude dependent and we do 
not find unambiguous identifications for the lower amplitude 
frequencies. The complete list of (l,m) identifications derived 
from the combination amplitudes are given in Table [4] 

Our analysis of EC14012-1446's combination amplitudes 
indicates that 9 of the dominant independent frequencies are 
consistent with spherical degree 1=1; in Table [4] we list the 
most likely {I, m) identifications for these modes. We would 
like to constrain the identifications of the remaining indepen- 
dent frequencies to improve our chances of success with non- 
linear convective light curve fitting. Using the combination 
analysis as a basis, we now turn to asteroseismology to pro- 
vide constraints on the identifications of the remaining fre- 
quencies. 

4.2. Asteroseismology 

4.2. 1 . Period Spacing 

As mentioned in the Introduction, we look for two clues 
to indicate (l,m) identifications of pulsations in white dwarfs. 
The first is the expectation that g-mode pulsations of a given I 
corresponding to successive radial overtones k will be approx- 
imately equally spaced in period, provided k is large enough 
( |Unno et al.| 1 989| >. For the simple example of a homogeneous 
: find 



star, we : 



P kl =kAH [l(l+l)]~ l/2 + C 



(1) 



where All is a uniform period spacing, F\i is the period fo r 
a given mode (k,l), and C is a constant (Tassoul et al. 1990). 
The mean period spacing All for a series of modes of a given 
I and consecutive k is an important asteroseismic measure of 
stellar mass and effective tem perature that is m ostly indepen- 
dent of internal composition (Unn o et al.|1989] >. 

A white dwarf has a host of available pulsation frequen- 
cies. For reasons that are not understood, in mos t cases only 
a subset are observed at any given time (Cor sico et al. 1 2002). 
This is true for EC14012-1446. The 19 independent frequen- 
cies detected during XCOV28 show no obvious evidence of 
equal period spacing so we are lacking consecutive radial 
overtones for any / value. A good strategy to identify the com- 
plete set of available modes is to combine results from mul- 
tiple seasons of observations. |Handler et al.| ([2008 ) present 
an analysis of EC14012-1446 observations spanning 2004- 
2007, during which the star exhibited different subsets of ex- 
cited modes. Figure 1 1 shows a schematic representation of 



excited modes for the combined observations spanning 2004- 
2008. Figure [12] focuses on the obvious groupings between 
800 and 500 s (1250 and 2000 /iHz). Since we are searching 
for equal period spacings, we present these figures using pe- 
riod as the x-axis. The large groupings at 768 s (1302 /iHz), 
721 s (1387 /xHz), 682 s (1466 /iHz), and 612 s (1633 //Hz) 
have a decreasing observed width with shorter period. This 
decrease translates into an equal width in frequency space of 
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Table 4 

Frequency (7,m) Identifications 



Frequency 


Period 


Combination 


Asteroseismology 


Final Nonlinear Fit 


fiHz 


(s) 


(l,m) 


(l,m) 


(/,m) 


935.3799 ±0.002 


1069.085 




(l.any) 


(2,-1) 


1021.139±0.002 


979.299 




(l.any) 


(2,2) 


1104.252 ±0.001 


905.591 




(l,any) 


(2,1) 


1155.925 ±0.002 


865.108 




(l,any) 


(2,2) 


1241.403 ±0.002 


805.540 




(l,any) 


(2,2) 


1418.369 ±0.002 


705.035 




(2,any) 


(2,1) 


1521.575 ±0.002 


657.214 




(1,1) 


1,1) 


1 548 146 + 001 


645 933 




(l,any) 


n o"i 


1624.015 ±0.003 


615.758 


(i,D 


(i i) 


(1,1) 


1633.907±0.001 


612.030 


(1,0) 


(1,0/1) 


(1,1) 


1774.989±0.100 


563.384 


(1,-1) 


(l,any) 


(1,1) 


1860.248 ±0.400 


537.563 




(2,any) 


(2,0) 


1883.555 ±0.003 


530.911 




(1,1) 


(1,1) 


1887.404±0.001 


529.828 


(1,0) 


(1,0) 


(1,0) 


1891.142±0.002 


528.781 


(1,-1) 


(1,-1) 


(1,-1) 


2304.745 ±0.001 


433.887 


(1,-1) 


(2,any) 


(1,0) 


2504.897 ±0.001 


399.218 


(1,1) 


(1,1) 


(1,0) 


2508.060 ±0.002 


398.715 


(1/2,1/0) 


(1,0) 


(1,-1) 


2856.155 ±0.002 


350.121 


(2,0/1) 


(2,any) 


(2,0) 



fa 24 [iHz. 

We calculated a simple average period for each group, 
which is given in Table [5J A statistical test for the presence 
of uniform period spacin g is provided by the Kolmogorov- 
Smirnov (KS) test ( [Winget et aLll99l) . The KS test calcu- 
lates the probability that an input list is randomly distributed. 
Any statistically nonrandom period spacing will therefore ap- 
pear as a minimum value in the output. In our case, we use 
the test to determine the probability that the list of average 
periods in Table [5] ar e from a unifor m distribution for a given 
period spacing All ( |Kawaler Jl988). Figure 13 shows the re- 
sults, finding a period spacing of 41 s, consistent w ith the 
expectations for 1=1 (Bischoff-Kim & Metcalfe 2011a). We 
find no significant period spacings at ~23 s as predicted for 
1=2 modes. 

Using 1633.907 /iHz as our reference point, we assigned 
relative radial overtone values to the frequencies listed in Ta- 
ble 5. We plot the expected locations of consecutive radial 
overtones for I = 1 in Figures 11 and 12 These points in the 



bottom of both figures are given with bars representing the 24 
/iHz "grouping" width. For periods above 900 s, the grouping 
widths for consecutive radial overtones overlap, illustrating 
the difficulty of assigning relative k values to modes with pe- 
riods longer than this value. Nonetheless, we find that most 
of EC14012-1446's pulsation frequencies can be identified as 
1=1. 

Unlike our simple example of a homogeneous star, white 
dwarfs are compositionally stratified, so the individual pe- 
riod spacings will not be uniform. We can retrieve detailed 
information about interior structure from the distribution of 
excited pulsation frequencies. A long standing problem with 
the asteroseismology of DA pulsating w hite dwarfs is the lack 
of objects w ith rich pulsation spectra (Bischoff- Kim & Met-| 
|calfe|j201 lb) . In this respect, EC14012-1446 immediately 
reveals its potential. Using our calculated average periods 
from Table I5J we show how individual period spacings AP 
differ as a function of period (and relative k value) in Fig- 
ure 14 We use "forward differencing", where AP is defined 
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Figure 11. Schematic representation of EC14012-1446's independent fre- 
quencies for data from 2004-2008. The bottom points indicate the expected 
locations of consecutive I = 1 modes with a mean period spacing of 41 s. 
The associated bars illustrate the "grouping width" of ±24 /xHz. Frequen- 
cies within this range should belong to the given k multiplet. Note that the 
.Y-axis is given in period, not frequency. 




Figure 12. A zoom-in of Figure [TT] to periods in the range 800-500 s. The 
bottom points give the expected locations of consecutive 1=1 modes with a 
mean period spacing of 41s. The period width denoted by the bars represent 
the 24 /iHz "grouping width" of the modes. Frequencies found within this 
range should belong to the given k multiplet. Note that the Ar-axis is given in 
period, not frequency. 
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Table 5 

Table of Average Frequencies (2004-2008) 



Period 


Frequency 


No. of peaks 


radial overtone 


(s) 








2856.164 


350.120 


1 




2738.076 


365.220 


1 


k -6 


2505.405 


399.137 


6 


k -5 


2304.854 


433.865 


2 




1887.505 


529.800 


10 


k -2 


1860.431 


537.510 


2 




1774.954 


563.395 


2 


fco-1 


1633.597 


612.146 


16 


k 


1548.155 


645.930 


2 


k + l 


1463.668 


683.215 


8 


k + 2 


1418.353 


705.040 


1 




1384.903 


722.072 


14 


k + 3 


1332.851 


750.249 


1 




1298.863 


769.904 


6 


k + 4 


1219.836 


819.782 


2 


k + 5 


1163.163 


859.710 


2 


k + 6 


1136.538 


879.865 


4 


k + l 


1 104.240 


905.600 


1 




1035.443 


965.770 


2 


k + 9 


964.748 


1036.540 


3 


ko + U 


935.375 


1069.090 


1 




821.254 


1217.650 


1 





Note. — This table gives the simple average frequen- 
cies for the groupings in Figure [TTJ and Figure [T2] The 
column "No. of peaks" gives the number of frequencies 
contained within each group. The radial overtone k is im- 
possible to determine observationally, so we give values 
relative to 1633 fiHz, defined to have k a . The relative 
k values denote a series of nearly consecutive modes of 
spherical index 1=1 with a period spacing of « 41 s. 

as AP = Pk-Pk+x- The filled points represent the periods be- 
tween pa 900-500 s that we are certain are consecutive radial 
overtones of 1=1. The open circles represent those periods 
above ss 900 s with ambiguous k identifications. The roughly 
cyclic behavior in Figure [T4"|is a sign of mode trapping. Mode 
trapping occurs in compositionally stratified stars when there 
is a resonance between a pulsation frequency and a surface 
layer. In theoretical models, a resonance occurs when a radial 
node coincides with a transition layer. In a DA white dwarf, 
transition zones occur at the boundary of the hydrogen and he- 
lium layers, the helium layer and the carbon/oxyg en core, and 
at poin ts in a possibly chemically stratified core (Montgomery 
2005b). The trapping cycle (the number of frequencies from 
minimum to minimum) is most sensitive to the location of the 
trapping layer, e.g., the base of the surface hydrogen layer. In 
addition, the trapping amplitude (the depth of the minima) is 
sensitive to the density gradient in the composition transition 
zone. In general, a larger gradient produces a larger trapping 
amplitude. Work is underway to deter mine EC14012- 1446's 
detailed ma ss and internal structure (Bischoff-K im & Met- 
|calfe|2011a| l. 



4.2.2. Multiplet Structure 

A second clue to constrain (l,m) identifications of pulsa- 
tions in white dwarfs is the presence of rotationally split mul- 
tiplets. The multiplet components have the same (k,l), and 
are further described by the azimuthal index m, which takes 
integer values between -/ and I. To first order, the frequency 
difference relative to the m = component of the multiplet is 
given by dviam = ~m £1(1 -Cm), where is the rotation period 
and Qj is a coefficient that depends on the pulsation eigen- 
functions evaluated in the nonrotating case. In the high-A: 



a -2 




40 

Period Spacing (s) 

Figure 13. Kolmogorov-Smirnov (KS) test applied to the average periods in 
Tablel5] A period spacing is defined as significant with a confidence level of 
(1-0X 100%. The results reveal an average period spacing of ss 41 s with 
a confidence level of 99.99%. 



asymptotic limit for g-modes, C# ~ 1 /t(£+ 1), although mod- 
els predict it to vary by w 10 percent over the range of ob- 
served periods in EC 1401 2. We adopt the convention that a 
positive value of m represents a progra de mode; retrograd e 
modes are denoted by negative m values (Winget et al. 1994). 
An additional important diagnostic is given by the ratio of 
splittings for /= ! and 1=2, where we expect 8vi = \ j 8v^ = 2 ~ 0.6 
for g-modes (Wing eTet al.|1991[ ). 

For EC140T2-1446 during XCOV26, we find one triplet 
(1887 ^Hz) and two doublets (1633 and 2504 /zHz) among 
our list of 19 frequencies (Table |2j. We begin by assum- 
ing that the 1887 /iHz triplet is a rotationally split multi- 
plet. Our combination analysis argues that the central com- 
ponent (1887.404 ^Hz, 529.828 s) is (1,0) and the 1891.142 
/iHz (528.781 s) component is (1,-1). The 1883.555 /xHz 
(530.911 s) component is low amplitude and is not identified 
as a parent of any detected combination frequencies. We can- 
not constrain its (l,m) identification via that method. Given 
(l,m) identifications for the first two modes, the 1883.555 




750 700 650 

Period (s) 



Figure 14. AP vs P for EC14012-1446, using the average periods from Ta- 
bled We use forward differencing (AP = Pj — Pfc+i ), where the radial over- 
tone/: increases to the left. The solid points represent those frequencies with 
1=1 identifications. The open points represent those where our identification 
is uncertain. 
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/iHz component should be (1, 1). The average multiplet split- 
ting of 3.8 /iHz implies a rotation period as sampled by these 
modes of 1.53 ±0.01 days. This is a perfectly reasonable 
result, as spectroscopic studies of white dwarfs reveal upper 
limits of v sin i = 10 k m s~ 2 , consistent with rotation periods 
of order days or longer (Berg er et al.|2005) . 

In the limit of uniform slow rotation, we expect additional 
/ = 1 multiplets to exhibit similar splitti ngs and 1= 2 modes to 
have splittings near 6.3 /iHz ( |Winget et al.|19 9lj). Recalling 
that every m component is not necessarily excited, we turn 
to the doublet at 1633 /iHz (612 s). Period spacing argues 
that this doublet is 1=1, and our combination analysis indi- 
cates that the 1633.907 /iHz (612.03 s) mode is (1,0). The 
1624.015 /iHz (615.758 s) component is not found as a par- 
ent of any detected combination frequencies. If 1633.907 
/iHz is indeed the central component of an /=1 triplet, then 
1624.015 /iHz is m = 1 and we do not detect the m = —l com- 
ponent. However, the splitting of 9.89 /iHz differs signifi- 
cantly from the expected value of w 3.8 /iHz found for the 
1887 /iHz triplet. A second possibility is that the combina- 
tion analysis is incorrect and we are missing the multiplet's 
central component. In this case, the 1633.907 /iHz mode is 
(1,-1) and the 1624.015 /iHz is (1, 1), resulting in an average 
splitting of 4.95 /iHz (9.89/2). Moving to the 2504 /iHz (399 
s) doublet, period spacing again argues that this multiplet is 
1=1. Our combination analysis indicates that the 2504.987 
/iHz (399.204 s) component is (1,1). Again, the 2508.060 
/iHz (398.715 s) component is not identified as a parent of 
any detected combination frequency. The frequency splitting 
between the two components is 3.16 /iHz and by compari- 
son with the splitting of the 1887 /iHz triplet, we can argue 
that these modes have consecutive m values, identifying the 
2508.060 /iHz component as (1,0). 

Multiplet splittings may also be used to eliminate the pos- 
sibility that EC14012-1446's frequencies represent a mix- 
ture of 1=1 and 1=2 modes. Assuming the 1887 /iHz 
triplet is 1=1, we do not find evidence for 1=2 splittings of 
5vi=2 = (3.8/0.6) =6.3 /iHz as predicted in the limit of slow, 
uniform rotation. Although both the period spacing and the 
combination analysis argues against it, to play devil's advo- 
cate we consider the possibility that the 1887 /iHz triplet is ac- 
tually an / = 2 quintuplet, since all m components are not nec- 
essarily excited to observable levels. In this case, then the ex- 
pected / = 1 multiplet splitting is 5vi = \ = 3.8 x 0.6 = 2.3 /iHz. 
We find no examples of multiplet splittings near this value. 
Finally, we must consider that Sv^^/Sv^ = 0.64. Based on 
multiplet structure alone, we could argue that the 1633.907 
/iHz frequency is I = 2, and 2504.897 /iHz is I = 1, but this 
is not supported by either the period spacings or the combi- 
nation analysis and leaves no clear explanation for the 1887 
/iHz triplet. 

Clearly, the multiplets found in EC14012-1446's XCOV26 
data set are not well explained based on the simple model 
of rotational splitting. Multiplet structure should be deter- 
mined by the star's rotation rate and structure. We expect 
this to remain unchanged over time. The classical example 
is PG1 159-035, which exhibit s triplets and quintuplets corre- 
sponding to I = 1 and 1 = 2 ( |Winget et al.|[T9 9T). However, 
complex multiplets are not unusual for white dwarf pulsators. 
In some instanc es, the ass umption of rigid rotation is clearly 
violated (Corsi ccTet al.|20lT) . In cooler pulsators with mod- 
erate amplitudes and well developed convection zones, ob- 
served multiplet structure can exhibit complicated behavior. 



I m=+l splitting 
I m=-l splitting 




1 2 

Relative radial overtone k 



Figure 15. Average multiplet structure associated with the large groupings 
at 770 s (1302 /iHz), 722 s (1387 /iHz), 612 s (1633 /iHz), and 399 s (1887 
/iHz). The m = +1 (prograde) splitting is always larger than the m = —1 split- 
ting, and the asymmetry increases with increasing k. 



For example, |Provencal et al.| ( |2009[ ) show changes in the 
DBV GD358's multiplet structure that cannot be e xplained 
with simple rotational splitting. For this same star, Winget 
et al. ( 1994) 1 show a dependence of multiplet splittings that 
is not explained by expected variations of Cu- Processes that 
may play a role in multiplet structure of cooler pulsators in- 
clude changing weak magnetic fields similar to the solar cy- 
cle, oblique pulsation, and differential rotation. Our XCOV26 
data on EC14012-1446 provides us with a single snapshot of 
this star's multiplet structure. We need more observations to 
understand their behavior. 

We turn again to the combined observations from 2004- 
2008. The four large groupings at 768 s (1302 /iHz), 721 s 
(1387_/iHz), 682 s (1466 /iHz), and 612 s (1633 /iHz) in Fig- 
all show indications of multiplet structure within their 
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groups. We extracted an average multiplet for each of the four 
large groupings in Figure 12 by calculating a simple average 
frequency value for the m =T, m = 0, and m = +1 components. 
The results are presented in Figure [15] In each case, the pro- 
grade (m = +1) mode splitting is larger than the m = -1 mode, 
and both the splittings and the asymmetries increase with in- 
creasing radial node k. 

The qualitative behavior of the low-A: and high-A: modes can 
be explained in terms of a very general model, g-modes are 
standing waves of buoyancy in a spherical cavity, and can be 
thought of as superpositions of traveling waves bouncing back 
and forth between an inner and outer turning point. In gen- 
eral, low-A: (shorter period) modes have deeper outer turning 
points while high-A: (longer period) modes have turning points 
closer to the stellar surface, meaning that these modes sample 
the outer regions of the star more than do low-k modes. Fig- 
ure [15] shows that the unknown process affecting EC14012- 
1446 7 s multiplet structure acts more dramatically on the high- 
k modes, arguing that the structural perturbation must be in 
the outer layers. A surface magnetic field and/or the convec- 
tion zone are obvious candidates. In addition, we should con- 
sider both radial and latitudinal differential rotation. Work is 
underway to improve our understanding of multiplet structure 
in pulsating white dwarfs (Dalessio 2012). 

Asteroseismology of EC 1 40 1 2-1 446 's distribution of ex- 
cited frequencies shows that the observed pulsations are dom- 
inated by a series of /=1 modes. If 1=2 modes are present, 
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Figure 16. Simultaneous fit of the periods of 19 modes (solid line) to the light curves from mcdof)80401, soar080403, and soar080405 (filled circles). Note the 
change in x axis for each plot. 

the modes are low amplitude and will not greatly affect the 
nonlinear light curve fitting process. For the 1887, 1663, and 
2504 /iHz multiplets, we are able to constrain the m identi- 
fications as well. The combination of these results with the 
combination amplitude analysis (Table HI gives us a strong 
foundation on which to proceed to nonlinear light curve fit- 
ting. 



5. NONLINEAR ANALYSIS 



Montgomery (2005a) and Montgomery et al. (2010b) give 
a detailed description of fitting observed non-sinusoidal light 
curves of white dwarfs to extract the time-averaged thermal 
response timescale of the convection zone. To summarize, 
the aspect of the convection zone sampled by the pulsations is 
the thermal response timescale, r c , which is directly related to 
the convection zone's mass and depth, and therefore its heat 
capacity. A DAV such as EC14012-1446 will experience lo- 
cal temperature excursions of ±250 K. In response, the local 
mass and depth of the convection zone will vary throughout a 
pulsation cycle. 

MLT predicts that r c should scale as 



T c PS To (T eS /T eS fl) 



(2) 



where T e s is the instantaneous effective temperature, r e ff,o 
is the equilibrium T e g, tq is the time-averaged convective 
timescale, and N describes the sensitivit y of r c to changes in 
r eff . For DAs, MLT pr edicts N^90 (Montgo meiy|2005a[|Wul 
2001 Brickhill 1992). The convection zone's mass, depth, 
and heat capacity are therefore extraordinarily temperature 
sensitive, and can vary by a factor of m 10 throughout a pulsa- 
tion cycle. This modulation of depth and heat capacity is the 



source of the large nonlinearities in white dwarf light curves. 
We focus on these observed nonlinearities to determine the 
time-averaged convective timescale to, the temperature sensi- 
tivity parameter N, and the inclination of the pulsation axis to 

the line of sight 6. 

Our a nalysis follows the approach of Montgomery et al. 
pOlOb) . EC14012-1446 is a multiperiodic pulsator, and be- 
cause of nonlinear effects, a pulse shape obtained by folding 
its light curve at a period of interest is not equivalent to a 
pulse shape obta ined in the absence of additional frequencies 
(Montgomery 2007). Therefore, we use the accurate frequen- 
cies and {l,m) identifications obtained for the 2008 WET cam- 
paign to calculate point by point nonlinear light curve fits to 
light curves obtained during XCOV26. High signal-to-noise 
light curves (S/N«1000) are vital, since we are interested in 
the nonlinear portion of the data, which is smaller than the 
linear component. We chose the SOAR 4m light curves, and 
the two longest runs from the McDonald 2.1m (selected runs 
are marked with "*" in Table [T] These 6 runs span three 
days, a timebase long enough to constrain the phases of the 
closely spaced frequencies, yet short enough to avoid possi- 
ble implications of amplitu de m odulation found in our fre- 
quency analysis (see Section 3.1 1. The dominant frequencies 



can be considered to be stable over a 3 day timescale. Using 
model atmosphere tables provided by D. Koester, we calcu- 
late the convers ion from bolometric to the observed bandpass 
as described in Montgomery et al. ( |2010b|) , assuming the fol- 
lowing parameters for this star: r e ffo= 11768 ± 23 K and log 
g=8.08 ± 0.008 in cgs units ( [Koester et al J2009| ). 
We began the fitting process by including the frequencies 
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in Table [4] with firm (7, to) identifications. We then included 
the additional frequencies with strong (/,to) constraints. It 
is prudent to point out that this fitting process is nonlinear. 
Adding additional frequencies based on criteria such as ampli- 
tude is not necessarily the best procedure. We experimented 
extensively by computing numerous fits encompassing a wide 
range of (I, to) identifications for the lower amplitude frequen- 
cies. The values of to for all fits range from 99-230 s, indicat- 
ing that tq is not strongly dependent on these identifications. 

Our investigation does reveal that the temperature parame- 
ter AT can be highly sensitive to the input to identification for 
large amplitude modes. Numerous attempts to fit EC14012- 
1446's light curve while assigning the 1633.907 /iHz mode 
an (7, to) identification of (1,0 ) as derived from our combina- 
tion analysis (see Section 4.1 1 yielded to=171 s, 8 = 28°, and 
N = 39. While the values of tq and 8 are reasonable values 
based on MLT, the value of N implies a tem perature sensitiv- 
ity that is far below predictions of MLT (Montgomery et al. 
2010b); this is not a physically relevant fit. We also experi- 
mented by assigning a spherical index 1=2 to this mode, re- 
sulting in the following fit parameters: to = 228 s, N = 26, 
8=15° . Again, the value of N obtained remains too low. |Ising| 
|& Koe ster (2001J warn that a perturbation analysis of combi- 
nation amplitudes dWu 2001) may have difficulties for pho- 
tometric variations of high amplitude, and this mode has the 
largest amplitude in the data set. In addition, the asteroseis- 
mology analysis of the 1633 /iHz multiplet structure c anno t 
solidly constrain the to value for this mode (see Section |4~2| i. 
Following this reasoning, we experimented by assigning the 
1633.907 /iHz mode an (l,m) value of (1,1). The resulting fit 
finds N = 85, in m uch better agreement with MLT predictions 
( |Wu|200 1 [ pontgomery|2005aj >. 

The final simultaneous nonlinear fit to the 6 high signal- 
to-noise light curves includes the frequencies and (l,m) iden- 
tifications given in Table 14] This fit produced the following 
parameters: r = 99.4± 12s, N = 85 ±6, 8 = 32.9±3.2°. Fig- 
ure [16] shows the ability to reproduce the essential features of 
the light curves. 

6. DISCUSSION 

Convective energy transport in stellar environme nts is typ- 
ically modeled u sing MLT. A particular version by Bohm & 
|CassinelE| ( |1971| ), denoted as ML2, includes re duced horizon- 
tal ene rgy lo ss relative to the formulation of Bohm-Vitense 
(ML1, |1958| l, increasing the overall convective efficiency. 
ML2 has been the standard convection model adopted for stel- 
lar atmosphere fits of white dwarfs for the past 20 y ears, with 
a = 0.6 the preferred value for the mixing length (Bergeron 
et al.| 1 995 ). Tremblay et al. (2010) recently re -calibrated the 



assumed convective efficiency for white dwarf models, using 
model spect ra incorporating an improved tr eatment of Stark 
broade ning ( |Tremblay & Bergeron 2009}. |Tremblay et al.| 
(|2010[ ) fit an improved set of Hubble Space Telescope and 
International Ultraviolet Explorer ultraviolet (UV) and near 
UV spectra, varying a until reaching agreement between the 
optical and UV temperatures. They find the best internal con- 
sistency between optical and UV effective temperatures and 
log g measurements using ML2 with a=0.8. This is a much 
more e fficient version of MLT than found by Bergero n et al.| 
( 1995 ), but is in closer agreement with that required by nona- 
diabatic mode ls to fit the observed blue e dge of the DA in- 
stability strip (Fontaine & Brassard 2008 i 



This convective 



parameterization is becoming the standard for DA model at- 
mospheres ( [Freytag et al.|2012||KllTc et al.|2012| >. 
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Figure 17. Comparison of EC14012-1446's derived convective parameters 
with values expected from ML2 convection. The addi ti onal labeled points 
are individu al objects taken f rom Montgomery (2005a I, Montgome ry et alT] 
1 2010a], and |Provencal et al.|(20ll) . The curves represent ML2 theoretical 
calculations of the thermal response time to for various values of the mixing 
length q (log g=8.0). 

One of our goals is to produce an empirical map of the time- 
averaged convective timescale to at the base of the convection 
zone as a function of T e g and logg for a population spanning 
the DAV instability strip (ss 1 1, 100-12,200 K) and compare 
this with MLT predictions of To. As an individual DA white 
dwarf cools through the instability strip, the base of its sur- 
face hydrogen convection zone will deepen, increasing its to- 
tal mass and the value of to. MLT makes specific predictions 
for the behavior of to given different choices of input param- 
eters. For instance, Figure 17 shows that ML2 with a=l.O 
predicts to values ranging from ~ 1 s at the blue edge of the 
instability s trip to ~ 6000 s at the red edge. 



In Figure 17 we also plot current determinations of Tq ver- 
sus T e g for five DAVs, including our solution for EC14012- 
1446; the derived stellar parameters are listed in Table [6] To 
ensure a uniform treatment of all the stars in our sample, we 
have used T e g determinations that do not includ e the recent 
updates to the line profile calculations by Tremblay & Berg- 



eron 



(2009 1; rather, the plotted effective temperatures and hor- 
izontal error b ars are based on earlier spectroscopic fits em- 
ploying ML2 ([Koester et al |2009j [Koester & Holberg|2001 



Koest er & Allard||2000| l. Overlaid on this fi gure are the to 
predictions of ML2 convection for various values of the mix- 
ing length parameter a. Our current results do indicate an 
increase in To (and hence an increase in depth and mass of 
the convection zone) with decreasing temperature, and are 
marginally consistent with a = 1.0. Decreasing the T e s error 
bars, either through higher signal-to-noise spectra and/or new 
model atmosphere fits, will provide more precise constraints 
on convection in these stars. 

Finally, we point out that convective light curve fitting ex- 
tracts the value of to at the base of the convection zone. There 
is no reason why the same value of a should describe both the 
photosphere and the deeper convective layers, so our results 
do not neces sari ly have to agree wit h the results of Tremblay 



et al.|(|2010[) or (Bergeron et al.|1995[). For instance, |Ludwig 



et al.| ( |T9"9"4"P use sophisticated 2-D hydrodynamic simulations 
to show that, while MLT is a reasonable approximation to pre- 
dict the rough photospheric temperature structure of DA white 
dwarfs, the deeper layers have a higher convective efficiency 
than predicted by MLT 



Whole Earth Telescope Observations of EC14012-1446 
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Table 6 

Stellar Parameters 



Star r eff (K) logg Period (s) t (s) P max 9 N 



G117-B15A 


12010± 180 


8.33 


215.2 


30 ± 10 


188 ± 63 






EC14012-1446 


11768 ± 23 


8.11 


612.3 


99 ± 12 


624 ± 75 


33 ± 3 


85 ± 6 


G29-38 


11690± 120 


8.11 


615.2 


187 ± 20 


1177± 126 


66± 3 


95 ± 8 


WDJ 1524-0030 


11660± 180 


8.06 


697.3 


163 ± 35 


1024± 220 


58 ±10 


95 ± 15 


GDI 54 


11270± 170 


7.95 


1185.9 


1169±200 


7345 ± 1257 


10± 10 


127 ±25 




with' 



: unique 



by l 



Our treatment of the nonlinearities in white dwarf light 
curves is based on the larger picture of how a surfac e convec- 
tion zo ne le ads to driving in these sta rs. In particular, Brickhill 
( |1991| ) and |Goldreich & Wu| ( |1999| > demonstrate that excita- 
tion of g-mode pulsations should occur when the convective 
driving exceeds the radiative damping, and that this condition 
is given by ujtq > 1 . In terms of period, this relation states that 
g-modes should be driven when P < P max , where P m . dx = 2htq 
and P is the period of a given mode. In Figure [18] we com- 
pare the dominant oscillation period in each of theDAVs we 
have fit with their theoretical value of P mdx , as calculated from 
each star's value of to- The agreement is good for tq < 100 s, 
but the values diverge for larger values of To. Given that tq 
is related to a star's T e g (see Figure 17 1, this says that the 
agreement is good from the blue edge to near the middle of 
the instability strip, but that from this point to the red edge 
another effect is operating that prevents the dominant mode 
periods from increasing as rapidly as P max . 

It is natural to assume that the cause of this behavior is 
related to whatever unknown phenomenon produces t he red 
edge of the instability strip. Wu & Goldreich d2001 ) point 



out that a deepening convection zone attenuates the flux vari 
ations entering at its base, so that for deep enough convection 
zones the surface amplitudes will be below detection limits. 
However, this is an amplitude effect and does not pr edict that 
the dominant periods should be different from P mix . |Brickhill 



u 
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Pmax 



Period (sec) 
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Figure 18. Comparison of the dominant observed pulsation period in a star 
(dot-dashed curve, red squares) with the theoretical value of P raax (dashed 
line, blue circles) for each DAV star with a measured value of to. Cooler 
stars are to the left, hotter stars to the right. While simple theory predicts that 
these curves should be very similar, we find a significant departure for cooler 
stars. 



( [T99l) , |GoIdreich & Wu| ( [T999] l, and |Wu & Goldreich| ( |200T] i 
also discuss the importance of turbulent damping in a narrow 
shear layer at the base of the convection zone; damping from 
this region may be more significant for cooler stars. Another 
possibility is that the maximum period for which complete 
surface reflection of g-modes can occur is much lower than 
expected; typical estimates place this number on the order 
of ~ 5000 s for / = 1 modes and ~ 3000 s for I = 2 modes 
(Hansen et al. 1985), although it is possible that improved 



treatment of the surface layers could alter these estimates. At 
any rate, this phenomenon is likely related to the deepening 
of the convection zone toward the red edge; we are exploring 
other phenomena that could lead to this behavior. 

7. CONCLUSIONS 

Until the advent of convective light curve fitting, our Sun 
was the only star with empirical constraints on its convec- 
tion zone depth. Determinations of the time averaged convec- 
tive timescale to and the temperature sensitivity parameter N 
can now potentially be obtaine d for any white dwa rf pulsator 
of moderate amplitude (Montgom ery et aT7||2010b[ >; approx- 
imately two-thirds of all white dwarf pulsators show signifi- 
cant nonlinearities in their light curves. Our ultimate goal is to 
map to as a function of T^q and logg for a population span- 
ning the instability strips of both the DAV and DBV white 
dwarfs. Such a map will provide important empirical con- 
straints on convection for white dwarfs and eventually other 
types of pulsating stars. 

We have taken the first steps in this direction with our in- 
vestigation of EC14012-1446. XCOV26 produced 308.3 h 
of data, and our analysis has identified 19 independent fre- 
quencies distributed in 14 multiplets. Combined with archival 
observations, we have identified a series of 1=1 modes with 
an average period spacing of 41 seconds. EC14012-1446 is 
now one of the few DAVs with over a dozen identified modes 
in its pulsation spectrum. The large number of modes means 
that asteroseismology can be used to provide constraints on 
its interior structure. 

Future goals focus on the determination of to and TV for ad- 
ditional white dwarf pulsators. An increased sample size will 
improve our empirical map of To as a function of r e ff and also 
allow us to further explore the observed behavior of P and 
^max as a function of To. Our goals also include the reduction 
of errors associated with spectroscopic temperatures. Convec- 
tive light curve fitting demands uniform treatment of convec- 
tion between the spectroscopic temperatures we choose and 
the temperatures we use to calculate the light curve fits. Since 
WD 1524-0030 was n ot included in the new temperatures of 
Gianninas et al. ( 201 l| l that incorporate recent updates to line 
profile calculations, we use a consistent set of older published 
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temperatures for the objects presented in this work. Incorpo- 
rating the new temperatures will be a first step towards reduc- 
ing the spectroscopic errors. We must recalculate our existing 
nonlinear light curve fits using these updated T^s a nd logg de- 
terminations. On average, the temperatures of Gianninas et al. 
( |201 1[ ) increase over older published values by « 500 K. Pre- 
liminary work indicates that the empirical value of to will not 
change significantly with such an increase. However, the the- 
oretical MLT predictions for tq will change; larger values of 
a will be needed to keep to and the convection zone depth the 
same for these higher values of T e ff. 
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